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Abstract 

We study a relation between droplet configurations in the bubbling AdS geometries 
and a complex matrix model that describes the dynamics of a class of chiral primary 
operators in dual N = 4 super Yang Mills (SYM). We show rigorously that a sin- 
glet holomorphic sector of the complex matrix model is equivalent to a holomorphic 
part of two-dimensional free fermions, and establish an exact correspondence between 
the singlet holomorphic sector of the complex matrix model and one-dimensional free 
fermions. Based on this correspondence, we find a relation of the singlet holomor- 
phic operators of the complex matrix model to the Wigner phase space distribution. 
By using this relation and the AdS / CFT duality, we give a further evidence that the 
droplets in the bubbling AdS geometries are identified with those in the phase space of 
the one-dimensional fermions. We also show that the above correspondence actually 
maps the operators of N = 4 SYM corresponding to the (dual) giant gravitons to the 
droplet configurations proposed in the literature. 
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1 Introduction 



Recently much attention has been paid to the half-BPS sector in the AdS/CFT correspon- 
dence [U El 0] due to the following two interesting observations. First, it was discussed 
in that the dynamics of the chiral primary operators of M = 4 super Yang Mills 
(SYM) on S 3 x R that correspond to Kaluza-Klein (KK) gravitons, giant gravitons [7J |H] 
and dual giant gravitons Ej in the supergravity side is described by a complex matrix 
model (complex matrix quantum mechanics), which is obtained by reducing M = 4 SYM to 
one dimension. It was pointed out in [6 that a holomorphic sector of the complex matrix 
model is related to one-dimensional ((l + l)-dimensional) fermions in the harmonic potential. 
In fact, the one-dimensional free fermions arise from a hermitian matrix model (hermitian 
matrix quantum mechanics) whose action is formally the same as that of the complex ma- 
trix model jTTj, and the number of the degrees of freedom of the hermitian matrix model 
is equal to that of the holomorphic sector of the complex matrix model. Second, as shown 
by Lin-Lunin-Maldacena (LLM) JUJ, the half-BPS solutions of type IIB supergravity that 
preserve R x 50(4) x SO (4) isometry are characterized by a single function satisfying a dif- 
ferential equation. The single function can be determined by giving an appropriate boundary 
condition in a two-dimensional subspace in ten-dimensional spacetime. Giving a boundary 
condition is equivalent to specifying shapes of droplets on the two-dimensional subspace. 
This is so-called bubbling AdS geometries (for further developments, see [12-32]). The afore- 
mentioned gravity states should belong to the above class of the half-BPS solutions. It is 
plausible due to the first observation mentioned above that these droplets are identified with 
those in the phase space of the one-dimensional free fermions. In fact, some evidences that 
support this remarkable identification have been reported jTUj H2 UH] • 

However, a precise map between the complex matrix model in the SYM side and the 
droplet configurations in the bubbling AdS geometries has not been given in the literature. 
The purpose of this paper is to find this map and give a further evidence of the above 
identification. Revisiting the analyses in j51 E] , we first show rigorously that the singlet 
holomorphic sector of the complex matrix model is equivalent to a holomorphic part of 
two-dimensional ((2 + l)-dimensional) free fermions in the spherically symmetric harmonic 
potential. It has been discussed [331 Ell that these two-dimensional free fermions are related 
to the one- dimensional free fermions. Another interesting work relevant for us is Ref.|13j. 
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The authors of JSI identified fluctuations of the droplet of the AdS*, x S 5 geometry that 
are responsible for the KK gravitons. Then, by representing the effective action for the KK 
gravitons in terms of these fluctuations, they gave a supporting evidence of the identification. 
Motivated by these works, we next establish an exact correspondence between the singlet 
holomorphic sector of the complex matrix model and the one- dimensional free fermions. 
Based on this correspondence, we find a relation of the singlet holomorphic operators of 
the complex matrix model to the Wigner phase space distribution of the one-dimensional 
fermions. Using this relation and the AdS / CFT duality, we give a further evidence that the 
droplets in the bubbling AdS geometries are identified with those in the phase space of the 
one-dimensional free fermions. We also show that the above correspondence actually maps 
the operators of M = 4 SYM corresponding to the (dual) giant gravitons to the droplet 
configurations of the supergravity side proposed in EH . Most of our results also hold 
beyond the classical limit, namely at finite N, so that we expect to gain a clue from them 
in extending the bubbling AdS business to finite N. 

The present paper is organized as follows. Section 2 is devoted to review. In section 

2.1, we refer to a class of chiral primary operators of A^ = 4 SYM on S 3 x R (multi trace 
operators) and their correlation functions that are analogues of the extremal correlators of 
M = 4 SYM on R A and we are concerned with throughout the paper. We also describe a 
complex matrix model that is used in calculation of these correlation functions. In section 

2.2, based on Refs. |35tlT2*l l31|. we summarize the results for the one-dimensional free fermions 
in the harmonic potential. In particular, we describe the classical phase space density and 
its quantum analogue, the Wigner phase space distribution. In section 2.3, we review the 
bubbling AdS geometries and its conjectured relation to the classical phase space of the one- 
dimensional free fermions. In section 3, we study the complex matrix model. We reduce the 
singlet holomorphic sector of the matrix model to the holomorphic part of two-dimensional 
free fermions in the spherically symmetric harmonic potential. We also develop the second 
quantization of these two-dimensional fermions. This gives an exact correspondence between 
the singlet holomorphic sector of the complex matrix model and the one-dimensional free 
fermions. In section 4, based on the correspondence, we find a relation of the singlet holomor- 
phic operators of the complex matrix model to the Wigner phase space distribution. Using 
this relation and the AdS/CFT duality, we give a further evidence that the droplets in the 
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bubbling AdS geometries are identified with those in the phase space of the one-dimensional 
fermions. In section 5, using the correspondence established in section 3, we construct in 
the Hilbert space of the one-dimensional fermions the states that correspond to the (dual) 
giant gravitons, and reproduce the droplet configurations of the (dual) giant gravitons in 
the supergravity side proposed in (HI CD] • Section 6 is devoted to summary and discussion. 
In appendix A, we present the results for harmonic oscillators in one and two dimensions, 
which we use in the main text. In appendix B, we give another derivation of the relation in 
section 4. 

2 Review 

2.1 Chiral primary operators and extremal correlators in J\f = 4 
SYM 

In this subsection, we review some facts about a BPS sector of M = 4 SYM. We are concerned 
with the chiral primary operators that take the form 

p 

{Ji,J2,-,J P } = Y[Tr(Z Ja ), (2.1) 

a=l 

where Z = ^{<j>i + ^2) with <pi and <p2 being two of six real scalars of M = 4 SYM. Note that 
these operators are singlet under U(N) transformations and holomorphic, namely depend 
only on Z. For simplicity, we denote a set {J 1: J%, ■ ■ ■ , J p } symbolically by {J}. We first 
consider correlation functions in M = 4 SYM on R 4 of the type 

((O^ 0) \y)yoV m \ Xl ) ■ ■ ■ OV (M) \x M )), y,x h ---, x M £ R\ (2.2) 

which are called the extremal correlators in the literature. The charge conservation requires 
that 

M 

j(o) = ^/M (2.3) 

r=l 

where 

PO Pr 

J(°) = ^J Q0 , = (2.4) 

clq=1 a r =l 
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There is a non-renormalization theorem on the extremal correlators telling that one can 
calculate them by using only the free part of the theory Next, we consider the chiral 
primary operators (j2.1j) in M = 4 SYM on S 3 x R. We restrict ourselves to the lowest 
Kaluza-Klein modes on S 3 because we are interested in half-BPS operators, so 0^ Jl,j2 ' '"' Jp ^(t) 
has only the t-dependence. (Linear combinations of) these operators can represent the KK 
gravitons, the giant gravitons and the dual giant gravitons jSJ Elj ■ We consider correlation 
functions of the type 



where t Q > t r (r = 1, • ■ • , M). These correlation functions are counterparts of the extremal 
correlators (|2.2jl on R 4 . It is natural to consider that the non-renormalization theorem also 
holds for the correlation functions (|2.5|) and one can calculate them using the free part of the 
theory. Reducing M = 4 SYM to the free part of Z(t) yields a matrix quantum mechanics 
defined by 



where Z(t) is an iV x N complex matrix and the path-integral measure is defined through 
norm in the configuration space of the matrix 



The potential term in the action arises from a coupling of the conformal matter to the 
curvature of S 3 and we have rescaled the field and the time appropriately. In this pa- 
per, we assume the non-renormalization theorem for the correlation functions ()2.5|) and 
concentrate on their calculation through the complex matrix model 1)2. 6j) . The condition 
t > t r (r = 1, • • • , M) in ()2.5|) is naturally understood from the following fact, which was 
pointed out in The extremal correlators ()2.2|) calculated through the free part of the 
theory agree with the correlation functions (J2.5j) calculated through the matrix model if all 
the propagators {Z i j(x)Zl l (y)) = SikSji/(x — y) 2 in ()2.2|) are replaced with those of the matrix 
model (Zij(t)Z^(to)) = 5 ifc 5 i; e i( *-* o) /2 that are valid for £ > t. 



((o^ C0) >(t )ro {jC1)} (ti)---o {J(M)} M), 



(2.5) 




(2.6) 



\dZ(t)\\ 2 = 2Ti(dZ(t)dZ\t)). 



(2.7) 
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2.2 One-dimensional fermions in the harmonic potential 

In this subsection, based on El EH] , we summarize the results for a one-dimensional 
system consisting of N non-interacting fermions in the harmonic potential. The classical 
one-body hamiltonian is given by h c i = ^p 2 + |g 2 . First, we describe a classical aspect of 
this system. A useful object in discussing the classical aspect is the classical phase space 
density u c i(p,q,t), which takes the values or 1. Regions in which u c i(p,q,t) = 1 are called 
'droplets'. The total area of droplets is equal to 2nhN (we put h = 1): 

< ^p-u d {p,q,t) = N. (2.8) 

It is convenient to introduce the polar coordinates (r, 0), 

q = rcos(p, p = rsin0. (2.9) 

Then, u c i(p,q.t) is expressed as 

u d (p,q,t) = 9(r(<l>)-r), (2.10) 

where f(0) is the boundary profile function, namely the equation r = f(0) represents a shape 
of a boundary of a droplet. Here we assume that there is only a single droplet and f(0) is 
a single-valued function. For example, the ground state is represented by a circular droplet 
f (0) = r with Tq = 2N. By using the classical equations of motion q = p, p = —q, it is 
easy to show that u c \ satisfies the equation of motion 

d d ( d d \ 

-^u d (p,q,t) = — Ud(rcos(f>,rsm(f>,t) = lq— -Pg^J u d{p,q,t). (2.11) 

Next, let us discuss a quantum aspect of the system. The quantum theory of the system 
is described by the second-quantized hamiltonian 

H = J dq^(q,t) {-\^ 2 + \q 2 ^j ^(q,t) (2.12) 

and a constraint 

' dq^(q,t)^(q,t) = N, (2.13) 



which turns out to fix the total number of particles to N. Here ip(q,t) is a fermionic field 
that satisfies the anti-commutation relation 

{^t)^\q\t)} = 5{q-q'). (2.14) 
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ip{q,t) obeys the equation of motion 

i ^) = W5 , f ),«] = (-I^ + I,^( 5 , t ). (2.15) 

It is seen from f!2.15|) and ()2.14|) that ip{q,t) is expanded in terms of creation- annihilation 
operators C\ and C n with an anti-commutation relation {C m , C^} = 5 mn as 

oo 

i/>(q,t) = Y,d n e- iE « t <p n (q), (2.16) 

n=0 

where E n = n + | and f n (q) is the normalized one-body wave function of the n-th ex- 
cited state, which is defined in (jA.4|) . The hamiltonian (j2.12|) and the constraint (j2.13p are 
expressed as 

oo 

ft = ^2 E n ClC n , 

n=0 

oo 

Y,C\C n = N. (2.17) 



n=0 

The left hand side of the second equation is the number operator, as promised. The ground 
state is 

\Q) = Cl f _ 1 C] v _ 2 ---Cl\0), (2.18) 

where |0) is the Fock vacuum defined by C n \0) =0 (n — 0, l,---). Excited states are 
obtained by replacing some of the C\ in with C^s with n > N — 1 in such a way that 
the total number of C^'s is preserved. 

An object that plays a crucial role in our analysis is the Wigner phase space distribution 
u(p, q, t), which is defined by 

u(p,q,t) = J dxe^(q + ^t)^(q-^t). (2.19) 

It follows from ()2.15|) and (|2.13J) that the Wigner phase space distribution satisfies the 
equation of motion 

d ( d d \ 

u{p,q,t)=\q—-p—\u{p,q,t) (2.20) 



dt \ dp dq 

and a constraint 

dpdq 



2tt 



u(p,q,t) = N. (2.21) 



There is another important constraint on the expectation value of u, u(p, q, t) = (T\u(p, q, t)\T), 
where |T) is an arbitrary state that satisfies the constraint (|2.13|) . u(p,q,t) satisfies a con- 
straint 

u*u(p,q,t) = u(p,q,t), (2.22) 
where * stands for the Moyal product: 

(A * B)(p, q) = A(p, q)e^^~^^B{p, q). (2.23) 
It is sometimes convenient to consider a Fourier transform of u(p,q,t): 

M= f d ^e-W-^u{p, q, t). (2.24) 



It immediately follows from (|2.20j) and (j2.21|) that u(a,/3,t) satisfies the equation of motion 

-gu{a,(3,t) = h—- a — ) f/(n. It). (2.2:)) 
and a constraint 

u(0,0,t)=N. (2.26) 

It is easy to show that u(a, j3, t) also satisfies the algebra 

[u{a, (3, t),u(a', /?', t)} = 2i sin -(a/3' - a'(3)u(a + a',f3 + (3', t). (2.27) 

In fact, u(a,(3,t) is characterized by (l2~2oT) . (gUBD and (l2~27l) . 

An important property of u(p, q, t) is that its expectation value u(p, q, t) is reduced to the 
classical phase space density u c i(p,q,t) in the classical limit. Actually, because the Moyal 
product becomes the ordinary product in the classical limit, the constraint ()2.22j) is reduced 
to an equation (u) 2 = u, whose solution is u = or 1. The equation of motion ()2.2()j) and 
the constraint (|2.21j) also hold for u, and these equation and constraint are equivalent to 
(B and (l2~Tll . 

2.3 Half-BPS geometries and free fermion droplets 

In this subsection, we review the LLM description of the AdS^ x S* 5 geometry and the 
half-BPS fluctuation modes in [TO! ITS]. 
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All half-BPS geometries of type IIB supergravity which preserve the R x SO (4) x 5*0(4) 
isometry are obtained by LLM ^U|. These half-BPS geometries are given by 



ds 2 = -h 



-2 



dt + V i dx i 



i=i 



<i?/ 2 + V] dx i dx i 



i=i 



+ ye G (i^ + ye~ G ^ , (2.28) 



/i = 2y cosh G , z = - tanh G , 

2/<9^ = eijdjz , ?/(<9iV,- - djVi) = e^dyZ , 

5 t = -^/ 2 e 2G , B t = -\y 2 e~ 2G , 

F = dB t A (oft + 1/) + + dB , F = dB t A (dt + V) + B t dV + dB , 



^2.29) 



dB 



1 3 - 



1 3 - 



2 



where a single function z(x , x , y) obeys the differential equation 



didiZ + yd y 



d y z 







^2.30) 



and characterizes the solutions. The function | — z on the y = plane, which we denote 
by w, takes the values or 1 and defines droplets as u c \ does. In the LLM description, the 
AdS$ x S 5 geometry is given by 



z(r,y;r ) 

Vr 



r 2 — Tq + y 2 



2\/(r 2 + rl + y 2 ) 2 - 4r 2 r 2 1 
1 



f2.31} 







V, = — 



r 2 + rg + 



a/ (r 2 + Tq + y 2 ) 2 — 4r 2 ? 



- 1 



Here r is identified with the AdS radius R through r = R 2 = \/2N and (r, 0) are the polar 
coordinates of the x\-x 2 plane. The droplet configuration (w(xi,X2)) is a circular droplet 
with the radius r in the (xi,x 2 ) plane (see Fig. |TJ). 

In [T3] a subclass of metric perturbation around the AdS$ x S 5 is discussed. Such per- 
turbation is obtained by considering small ripples of the circular droplet. The boundary of 
the perturbed droplet in the polar coordinates is given by 



r (0) = r o + Sr((j)). 
Sr((j)) is expanded in Fourier modes as 

5r(0) = E a ^ < : 



[2.32) 



(2.33) 





Figure 1: Circular droplet corresponding Figure 2: Fluctuation around the circular 

to AdS 5 x S 5 . droplet. 

(see Fig. |2J. Here the zero mode is absent due to the area preserving deformation (to the 
first order in 5r). Then one can express the resultant metric perturbation in terms of a n . The 
result shows that the gravity modes described by a n are the KK gravitons that, according 
to the standard AdS/CFT dictionary, correspond to TrZ n for n > and Tr(Z^ n ') for n < 
in dual M = 4 SYM. Moreover, one can read off the commutation relation among a n 's from 
the effective action that gives the equations of motion for this metric perturbation: 



This commutation relation agrees with the one that is obtained when the droplet is regarded 
as the droplet on the classical phase space of one-dimensional fermions in the harmonic 
potential, with the identification x\ = q, x<i = P- Thus we are allowed to identify droplets 
in the x\ — x^ plane with those in the (p, q) phase space of the one-dimensional fermions. 

3 Complex matrix model 

3.1 Complex matrix model and free fermions in two dimensions 

In this subsection, we study the complex matrix model introduced in section 2.2. In particu- 
lar, we show that if we restrict ourselves to the correlation functions f)2.5j) . we can regard the 
model as a holomorphic part of a two-dimensional system of free fermions in the spherically 
symmetric harmonic potential. The action in ()2.(jj) is evaluated as 



s = / dm{z(t)z\t) - z{t)z\t)) = / d* £(z(f) y z(f); - zitu^zit)^, (3.1) 



[a m , a n ) ~ rnS,, 



m+n ■ 



(2.34) 
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while the norm (|2.7|) as 

\\dZ{t) || 2 = 2Ti(dZ(t)dZ\t)) = 2 d z {t)ijdZ{t)*ij, (3-2) 

from which we find an explicit form of the path-integral measure: 

\dZ(t)dZ\t)} = [2 N2 WdReZifyijdlmZit)^}. (3.3) 

hi 

By comparing ()3.1|) and (j3.3|) with ()A.15|) and (jA.16|) in appendix A, we see that the system is 
nothing but a set of A^ 2 independent two-dimensional harmonic oscillators that are spherically 
symmetric. 

We can quantize the system canonically as we do for the one-body system in appendix 
A. In what follows, we make use of the results in appendix A. The quantum hamiltonian is 

and the normalized ground state wave function is 

Xo = 4 re -^ t ) = 4 ? e-^^. (3.5) 

7T 2 n 2 

The measure for the inner product between the wave functions is 

J Y[ dZijdZ-j = 2 N2 f]J dReZijdlmZij. (3.6) 

ij ij 

We see from the results in appendix A that wave functions that take the form 

X {Ji '---' Jk) = (j[Tr(Z J ^ Xo (3.7) 

are eigenf unctions of H with eigenvalues A^ 2 + J2b=i Jb- This observation enables us to 
calculate the correlation functions ()2.5|) as follows: 

((^ (O) >(t O ))^ {j(1)} (tl)---0 {j(M)} (tM)) 

= (xo|e^°(C {J(0)} )*e~^ ( * 0_ ' l) C {J(1)} e _ ^ ( ' 1_ * 2) • ■ ■ e ~ i{l{tM '^ tM) {j{M)} e~ itltM |xo) 

= e ^L^(tr-t0) /]J^Z*.(0{ J(O ^^ (3.8) 

J ij 
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Here the result is independent of the order of ti, • ■ • , £m although we assumed in the above 
calculation t > t\ > ■ ■ • > tu- Hence it is necessary and sufficient for calculating the 
correlation functions ()2.5|) to examine the wave functions ()3.7|) . which are holomorphic except 
for the factor e~^- Ji ^ ZiiZ ^ and singlet under U(N) transformations. 

Let us see that as far as the wave functions (|3.7|) are concerned, we can reduce the 
dynamical degrees of freedom of the system to N eigenvalues of Z. We express Z in terms 
of a unitary matrix and a triangle complex matrix as 

Z = UTU\ (3.9) 

where UW = 1, Ty = for i > j and the Tu (i = 1, • • • , N) are the eigenvalues of Z |3T] . 
Note that the number of real parameters in Z is 2N 2 while those in U and T are iV 2 and 
iV(iV + 1), respectively. This reflects the fact that Z does not determine U and T uniquely: 
if V is a unitary diagonal matrix, (UV)^Z{UV) = V^TV and V'TV is also a triangle matrix. 
This redundancy enables us to impose N constraints on variation of U. Variation of Z is 

dZ = U(dT + i(dHT -TdH))U\ (3.10) 

where dH = —iU^dll and dH is a hermitian matrix. We impose the constraint dHa — (i — 
1, ■ • • , N). Because ||<iZ|| 2 = \\dT + i(dHT — TdH)\\ 2 , we can make a change of variables 
from Zij to Hij (i > j) and Ty (i < j). The jacobian for this change of variables is |A(z)| 2 
[3*7] . where 

A(z) = Y[(Zi-Zj), Zi = T u (i = l,...,iV). (3.11) 

i<j 

Hence we obtain 

I \\dZ, ] dZ] i = I HdH.jdH^lldTudT^lldz^JAiz)^ (3.12) 

i,j i>j k<l m 

This leads us to redefine a wave function x as 

XF = A(z) X (3.13) 

and the hamiltonian H as 

H F = A(z)H-^ (3.14) 
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Thus the system is equivalently described by the hamiltonian Hp and the wave functions 
Xf with the inner product 



(3.15) 



i>j 



k<l 



In particular, if x is an eigenstate of H with an eigenvalue E, xf is an eigenstate of Hp with 
the same eigenvalue E. Xf 1 ' "' Jk ^ is expressed in terms of Zj and (i < j) as 



A-h,-,J K ) 

Xf 



Xo 



\a=l i a / 



Xo, 



(3.16) 



N 2 
7T 2 



It should be noted that this eigenstate is independent of Hij and iif? due to the singlet 
nature of the wave function (|3.7|) . while this eigenstate is independent of T^-, T*j, and z* 
except for xo due to the holomorphy and the singlet nature of the wave function (|3.7|) . We 
see from this expression that xP' ' Jk ^ is equal to a certain linear combination of 



$^(21,4) ^(^2,^; 

$fe(*i>*i) 



$i N (z N ,z 



N) 



j<k 



1 K 
with = AT(iV-l) + ^J 6 



(3.17) 



6=1 



where &i(z, z*) is a 'holomorphic' wave function or a wave function of the lowest Landau level, 
which is defined in (jA.13j) . Namely, x^ 1 ' ' "' is a wave function of an energy eigenstate of a 
two-dimensional system consisting of iV fermions and |iV(iV — 1) bosons which are governed 
by the spherically symmetric harmonic potential and do not interact each other. In this 
picture the coordinates of the fermions are represented by (zj, z*), while the coordinates of the 
bosons are represented by (I}fc, IJJ. Here, the fermions are in the "holomorphic" states, and 
the bosons are in the ground state. The total energy of the fermions is |iV(iV" + l) + J2b=i 
where ^N(N + 1) corresponds to the energy of the ground state of fermions, and the total 
energy of the boson is ^N(N — 1). The sum of these two total energies is equal to the 
eigenvalue of Hp with respect to Xf ' ■> N 2 + Yl^=i Jb, so that we can identify the system 
with the system of these free fermions and bosons, as far as we concentrate on the wave 
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functions ' ■ Because the bosons are always in the ground state, their contribution to 
the eigenvalues of Hp is always ^N(N-l), and in the inner product (x F ' Kl \Xf 1 ' K2 ) 
the integral over Hij and 7y gives an overall constant factor. Therefore we are allowed to 
consider only the fermionic part of Xf' ' anc ^ re pl ace the measure in the inner product 
fl3~T5jl with jUidzidz*. In this way, we have reduced the singlet holomorphic sector of 
the complex matrix model to the holomorphic part of two-dimensional free fermions in the 
harmonic potential. This holomorphic nature enables us to relate these two-dimensional 
fermions to the one-dimensional ones, as we will see in the next subsection. 

Finally we make a comment on the singlet non-holomorphic sector which is non-BPS. The 
hamiltonian ()3.4|) describes a set of N 2 non-interacting harmonic oscillators represented by 
(Zij,Z*j). After change of variables from (Zij,Z*j) to (zi, z* ,Tij,T*j, H^, H*j), the hamilto- 
nian is described by the dynamical degrees of freedom represented by (zi, z*,Tij, T*-, H^, H*-) 
that are complicatedly interacting. Nevertheless, as we showed above, the singlet holomor- 
phic sector is reduced to the non- interacting fermions and bosons represented by (zi, z*) and 
(Tij,T*j) respectively. If we consider the singlet non-holomorphic sector, such simplification 
does not occur. We need to treat the above hamiltonian with the complicated interactions. 

3.2 Second quantization of two-dimensional fermions 

In this subsection, we develop a second quantized theory of the holomorphic part of the two- 
dimensional free fermions in the previous subsection and establish an exact correspondence 
between the singlet holomorphic sector of the complex matrix model and the one-dimensional 
free fermions. The observation in the previous subsection leads us to introduce a fermion 
field 

oo 

*(z,z*,t) = ^dte-^Qfaz*), (3.18) 
1=0 

where E\ — I + 1, and Ci and C\ satisfy the anti-commutation relation {C^C^} = 5i m . 
The C\ in (|3.18J) will be identified with C\ in ()2.16|) below. This fermion field satisfies an 
anti-commutation relation 

oo 

z*,t), *V, z'*,t)} = J2 ^ Z *)^' = K & z *~> z '' ( 3 - 19 ) 

1=0 
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K behaves as the delta function with respect to functions spanned by the &i(z, z*), namely 
holomorphic functions of z times e~ zz * . ty(z,z*,t) also satisfies a constraint that the total 
number of fermions be N, 

/oo 
dzdz*^(z,z*,t)^(z,z*,t) = J^CjCi = N. (3.20) 
1=0 

The ground state \ fl) is again given by 

|^) =C7^ 1 C'J v „ 2 ---C' 1 t C' t |0), (3.21) 

where |0) is the Fock vacuum defined by Cj|0) = 0. 
We define the operators §j(t) by 

§j(t) = J dzdz*z J U(z,z*,t), U(z,z*,t) = &(z,z*,t)V(z,z*,t), (3.22) 

where U(z,z*,t) is interpreted as the density operator restricted to the "holomorphic" sec- 
tor. It is easy to verify that one can calculate the correlation functions (|2.5jl by making a 
correspondence in the second line of (|3.8)l : 

Tr(Z(t) J ) <- sj(t) 

IXo) <- (3.23) 
The explicit form of ij(t) is easily calculated as 



*At) = ^ Jt Y,y^^ d lA- (3-24) 

Note that because [sj(t), sj>{t)] = 0, there is no ordering ambiguity in translating the multi- 
trace operators into products of sj's. Thus, if we identify C\ and C\ in §j with those 
of the one-dimensional fermions, we establish an exact correspondence between the singlet 
holomorphic sector of the complex matrix model and the one- dimensional free fermions in 
the harmonic potential. 

As an example, we give the calculation of (Tr(Z''"(to) 4 )Tr(Z(t 1 ) 2 )Tr(Z(t2) 2 )) with t > 
t\, t 2 , based on the formalism we developed above. 

(Tr^t^) 4 )^^^) 2 )^^^) 2 )) = (n\sl(t )s 2 (ti)§2(t2M 

= ^e 2iitl - to)+2i{t2 - to) (n\PQR\n), (3.25) 
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where 



Q 



R 



y/(N + 3)(N + 2)(N + l)NC j N ^C N+3 + y/(N + 2)(N + 1)N(N - l)C ] N „ 2 C N+2 
+ y/(N + PPV - l)(iV - 2)C*j v _ 3 CV + i + ^(iV - l)(iV - 2)(iV - 3)Ci_ 4 CW 
^(N-l){N-2)& N _ 1 C N - 3 + y/(N-2)(N-3)& N _ 2 C N -4 
+ y/(N + 3)(N + 2)Cl +3 C N+1 + y/{N + 2){N + l)& N+2 C N , 
vW+WcJr+iCW-i + V 7 ^ - 1)^^-2- (3.26) 



From ()3.25|) . we obtain the final result 

16iV 3 + 8iV 
2^ ' 



,2i(t 1 -t )+2i(t 2 -to) 



(3.27) 



which we can also obtain by applying the 'Wick theorem' with respect to Zij and Z*j to the 
last line of (J3.8|) . Note that the final result contains the non-planar contribution (See Figs. 

EE}. 1 



O.O 




16iV 3 



8A r 



Figure 3: planar contribution. 



Figure 4: nonplanar contribution. 



4 Singlet holomorphic operators 

and the Wigner phase space distribution 

In this section, based on the exact correspondence between the singlet holomorphic sector of 
the complex matrix model and the one-dimensional fermions, which was established in the 
previous section, we find an exact relation of the singlet holomorphic operators of the complex 
matrix model to the Wigner phase space distribution of the one- dimensional fermions. The 
relation gives via the AdS/CFT duality an evidence that the droplets in the bubbling AdS 
1 H3.27|l is valid for U(N) gauge group. For SU(N), the subleading term is modified |38|. 
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geometries are identified with those in the phase space of the one-dimensional free fermions 
in the harmonic potential. 

In section 2.3, we saw a n (n > 0) in ()2.33|) . which is the Fourier mode of the small 
fluctuation around the circular droplet (FigHJ), yields the KK graviton that corresponds to 
Tr(Z(t) n ) in dual Af = 4 SYM according to the standard AdS/CFT dictionary. In this case, 
w in section 2.3 that specifies the droplet is given by 

w(xi,x 2 ) = 9(r + 5r (ft)— r) 

= 6(r -r) + 5(r -r)5r((j)) + O((5r) 2 ), (4.1) 

where x\ = rcos0, x 2 = rsin</>. Thus, we obtain from ()2.33|) a relation 

1 dr \ d<f)rg n {r)e- in(i 'w{x 1 ,X2). (4.2) 



2nrog n (r ) 

At this stage, g n {r) is arbitrary. In the previous section, we constructed Sj(t) which corre- 
sponds to Tr(Z(t) J ) quantum mechanically. The quantum analogue of the classical phase 
space density u c i is u. Then, if the droplet specified by w(xi,x 2 ) are identified with that in 
the phase space of the free fermions through w(xi,x 2 ) = u c i(p,q) with x± = q and x 2 = p, 
we expect the following identity to hold quantum mechanically: 

s,j{t)=Aj [dr [ d<t>rgj{r)e- iJ(l> u(p,q,t). (4.3) 



Here an overall coefficient Aj and a function gj(r) will be determined soon. The existence 
of the relation ()4.3|) serves as an evidence that the droplets in the bubbling AdS geometries 
are identified with those in the phase space of the one-dimensional free fermions. This is one 
of our results in this paper. 

In what follows, we will see that the relation ()4.3j) indeed holds. We first evaluate it(p, q, t) 
explicitly. From (|2.16|) and (|2.19|) . we obtain 

oo „ 

fi(p,9,t) = E ClC n e^-^l dxe*»^(g + |K(9-|). (4.4) 

m,n=0 

In order to calculate the integral in (|4.4jl , we deform the contour as x — > 2x + 2ip and make 
use of the following formula: 

J dxe~ x2 H m {x + v)H n {x + w) = 2 n ^m\w n - m L^ m {-2vw) for m < n, (4.5) 
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where L™ m is the Laguerre polynomial defined by 



m=0 

The final result is 



w = EH)i" + ! j («) 



n — m I ml 



(p, ? ,t) = V (7t(7 n e i(m - n) S/^2^ +1 (-l) m e 4(m - n ^e- r r n - m L"- m (2r 2 ) 
z — ' V n! 

m<n 

+ ^ ClC n e^ m - n)t J 1^2^ + \-l) n e i{m - n ^e- r2 r m - n L™- n {2r 2 ). (4.7) 



TO! 

m>n 



By comparing (IT241 with (Oil and (I4~7f) . we find that (jOjl holds if Aj = and gj(r) 
r J . Namely, we find 



sj(t) = J dzdz*z J U(z,z*,t) 

dr I d<f)r J+1 e- iJ(t, u{p,q,t) 



1 



2tt2 j _ 

From ()4.8|) . we naively expect that a relation U(z,z*,t) = -u(p,q,t) holds with the identi- 
fication z = ~(q — ip). Unfortunately, this is not the case. It is, however, known jSl] that 
there is a simple relation between u(a,f3,t) and a Fourier transform of U. In appendix B, 
we give another derivation of ()4.8|) based on this relation. 

5 Mapping (dual) giant gravitons to droplets 

In this section, using the correspondence between the singlet holomorphic sector of the 
complex matrix model and the one-dimensional free fermions, which was established in 
section 3, we construct in the Hilbert space of the one- dimensional fermions the states 
corresponding to the (dual) giant gravitons. We will see that the classical phase space density 
that is the classical limit of the expectation value of the Wigner phase space distribution 
with respect to these states actually agree with the corresponding droplet configurations in 
the supergravity side proposed in [TU] . 

We first consider a sequence of non-negative integers (A) = (Ai, A2, • • • , Ajy) with Ai > 
A 2 > ••• > Aat > and £) i=1 Aj = n. This sequence defines a Young tableau with the 
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number of boxes in the z-th row being Aj and the total number of boxes being n, and 
specifies an irreducible representation of GL(N, C) as well as an irreducible representation 
of the symmetry group S n . As in [5], we introduce an operator which is a linear combination 
of the chiral primary operators 1)2.10 : 



X(\){Z) = ^Y. X(A)(<T)tr(<rZ), 



(5-12 



cres„ 



where X(\){Z) an d X(A)( cr ) are characters of the representations specified by (A), and tr(crZ) 
is defined by 

N 

tx(aZ) — 4v(i)4»„(2) ' ' ' (5- 2 ) 



H,i2i—|in=l 



As in [5 , we make use of Weyl's character formula. Weyl's character formula tells that 1)5.1)) 
is equal to the Schur polynomial: 



X(\){Z) 



det M zf-^ 
A(z) 



(5.3) 



Hence a fermionic wave function that is obtained by acting the operator (j5.1j) on the ground 
state is given by 

$N-l+\i( z li z *) ^JV-l+Ai ( z 2, z *2> ■'' ^JV-l+Ai (^AT, ^;v) 
$V-2+A 2 ( 2; l ) ^l) : 



Xf 



x]l$ (T jk ,T* k ). (5.4) 



$ Aflr (zi,zJ) •■■ •■■ <$> Xn (zn,z* n ) 

The correspondence established in section 3 maps this wave function to a state in the Hilbert 
space of the one-dimensional free fermions 

Av-l+AiAv-2+A 2 ' ' ' ^I+Aat-i^AatIO)' (5-5) 

This state can also be obtained by rewriting (j5.1j) in terms of sj through (J3.23)) and acting 
the resultant operator on the ground state \ 

Next, let us construct in the Hilbert space of the one-dimensional fermions the states 
that represent the (dual) giant gravitons. In M = 4 SYM side, the giant graviton with the 
angular momentum L, where L is compatible with N, corresponds to the operator ()5.1j) with 
(A) = (Ai = 1, A 2 = 1, • • • , \l = 1, Al+i = 0, • • • , Aat = 0) [HI 15] ■ This (A) corresponds to the 



18 



Young tableau that consists of a single column with L boxes. Note that there is a restriction 
L < N and the giant graviton with L = N is called the maximal giant graviton. We see from 
()5.5|) that the state that is obtained by acting this operator on the ground state is mapped 
to a state in the Hilbert space of the one- dimensional fermions 

\GG;L) = & N C N - L \n). (5.6) 

The phase space density that is the classical limit of (GG; L\u(p, q, t)\GG; L) obviously speci- 
fies the circular droplet with a small circular defect inside the circle (see Fig. |5J). This droplet 
configuration is indeed proposed in [SJ [TU] as that of the giant graviton in the bubbling AdS 
geometries. 




The dual giant graviton with angular momentum L, where L is compatible with or much 
larger than N, corresponds to the operator (j5.1j) with (A) = (L, 0, • • • , 0) 0. This (A) 
corresponds to the Young tableau that consists of a single row with L boxes. We again see 
from (|5.5|) that the state that is obtained by acting this operator on the ground state is 
mapped to a state in the Hilbert space of the one-dimensional fermions 

\DGG ] L) = C ] N _ 1+L C N _ l \Q). (5.7) 

This is also consistent with the proposal in [SJ ED] that the dual giant gravitons are given 
by the circular droplet and a small droplet outside the circular droplet in the bubbling AdS 
geometries (see Fig. [HJ). 
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6 Summary and discussion 

In this paper, we study a relation between the droplet configurations in the bubbling AdS 
geometries and the singlet holomorphic sector of the complex matrix model that can describe 
a class of the chiral primary operators of M = 4 SYM. We first showed rigorously that the 
singlet holomorphic sector of the complex matrix model is equivalent to the holomorphic part 
of the two-dimensional fermions. We developed the second quantization of these fermions and 
established the exact correspondence between the singlet holomorphic sector of the complex 
matrix model and the one- dimensional free fermions in the harmonic potential. Next, based 
on this correspondence, we find a relation between the operators Tr(Z(t) J ) of the complex 
matrix model and the Wigner phase space distribution. This gives via the AdS/CFT duality 
a further evidence that the droplets in the bubbling AdS geometries are identified with those 
in the phase space of the one-dimensional free fermions. Finally, we also constructed the 
states corresponding to the (dual) giant gravitons in the Hilbert space of the one-dimensional 
fermions. We obtained the droplet configurations for these states, which are consistent with 
the proposals in [HI ED] • Our main results are summarized in (jri.2Hj) , (jH.24j) , ()4.8j) , (|5.fi|) and 
(|5.7j) . These relations also hold at finite N, so that we expect to gain a clue from them in 
extending the bubbling AdS business to finite N. 

We make a comment on a hermitian matrix model whose action is formally the same as 
the complex matrix model: 



where the path-integral measure is defined by the norm ||<i0(t)|| 2 = Tr((c?0(t)) 2 ). As is well- 
known jllj . as far as the ?7(iV)-invariant sector is concerned, the model can be reduced to a 
system of N eigenvalues and these eigenvalues behaves as one-dimensional non-interacting 
fermions in the harmonic potential, which were discussed in section 2.2. Then, a U(N)- 
invariant operator Tr(0(t) J ) which is a counterpart of Tr(Z(t) J ) in the complex matrix 
model is translated to a operator in the second quantized theory of the fermions 



We cannot, however, replace sj with vj in the calculation of the correlation functions (|2.5|l . 
This yields extra terms that are formally obtained by 'contracting' Z^- and Z^. This is a 
manifestation of the difference between the complex and hermitian matrix models. 




(6.1) 




(6.2) 
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There are some directions as extension of the present work. First, it has been discussed 
that the pp-wave limit of the bubbling AdS space is described by relativistic fermions in 
(1 + 1) dimensions PHIED!- It is important to find a relation between these fermions and the 
complex matrix model. Topology changes in the bubbling geometry are also an interesting 
subject We hope our findings in this paper to give some insight to this subject. 

Acknowledgements 

We would like to thank Y. Susaki, T. Takayanagi and K. Yoshida for discussions. The work 
of Y.T. is supported in part by The 21st Century COE Program "Towards a New Basic 
Science; Depth and Synthesis." The work of A.T. is supported in part by Grant-in-Aid for 
Scientific Research (No. 16740144) from the Ministry of Education, Culture, Sports, Science 
and Technology. 

Appendix A: Harmonic oscillators in one and two di- 
mensions 

In this appendix, we summarize the results for the harmonic oscillators in one and two 
dimensions, which we use throughout this paper. The harmonic oscillator in one dimension 
is defined by a hamiltonian, 

t 1 ^ 1 2 ~t~ 1 ,A -,N 

h = -2W + 2 q = aW 2' (A ' 1} 
where a) and a are of course the creation and annihilation operators: 

1 / d , x 1 . d . . . . 

= vi (? W (A ' 2) 

which satisfy the commutation relation [a,tv] = 1. The ground state wave function <po(q) is 
characterized by a relation a<j) (q) = and the normalized one takes the form 

Md) = 4e~^ 2 . (A.3) 

The normalized wave function of the n-th excited state, whose energy eigenvalue is n + |, is 
explicitly given by 

<p n {q) = ^=(« t )>o(g) = * H n (q)e-h\ (A.4) 

vn! V7T22 n n! 
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where H n is the Hermite polynomial. 

We generalize the above results to the two-dimensional spherically symmetric harmonic 
oscillator, which is defined by a hamiltonian, 

i 1 92 1 92 1 2 1 2 n 

* = ~2^f "2^f + 2* + 2* 2 = ° iai + 4 ° 2 + 1; (A ' 5) 

where 

a,= 7f (?i+ l; ) ' 81 = - (t=1,2) - (A - 6) 

These creation and annihilation operators satisfy the commutation relations 

[a h a]] = Sij, [at, a,j] = [a], aj] = 0. (A. 7) 

The normalized wave function for the (m, n) excited state, whose energy eigenvalue is m + 
n + 1, is given by <p m , n (qi, q 2 ) = y?m (<?i)y?n(<?2 )• 

It is crucial for our analysis in the main text to introduce the complex coordinate: 

z=-j=( qi + iq 2 ). (A.8) 

In this variable, the hamiltonian takes the form 

d 2 

h = + zz * = 6 i fi i + 4c2 + 1, (A.9) 

ozoz* 

where 

1 d 

ci = -^(ai + ^2) = -^=(2 + cl = ~^{a\ - ia\) = -j=(z* - — ), 

+ i4) = 7! (2 -^ ) ' (A - 10) 

which satisfy again the commutation relations [q, cj] = 5^, [q, c 3 -] = [cj, cj] = 0. By acting 
the creation operators c\ and c| successively on the ground state, one can construct energy 
eigenstates, which form a normalized orthonormal basis: 

<M*>0 = i$W4)W- (A.11) 

yk\l\ 

The ground state, whose energy eigenvalue is one, takes the form in the complex variable 

S ,o(*,**) = 4= e ""*- ( A - 12 ) 
V 71 " 

22 



1 r 




' m 9 \ 


21 = 


1 
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1 r 

T2 {ai 


~ lh2) = T2 


< * 9 \ 

- + &)• 


4 = 


1 

7! 



The excited state &k,i{z,z*) with the energy eigenvalue k + I + 1 is expressed as a linear 
combination of ip m ,n with m + n — k + I. In particular, $ ,z is holomorphic except for the 
factor e~ zz ": 

$ Q j(z 1 z*) = ^Je-**\ (A.13) 

We denote $ j(z,z*) simply by $>i(z,z*). $i(z,z*) is a wave function of the lowest Landau 
level in the context of the quantum Hall effect. The measure used for the inner product 
between the wave functions is 

J dqidq 2 = 2 J dRezdlmz* = J dzdz*. (A. 14) 

The system can be quantized equivalently by path-integral. The partition function is given 
by 

Z = J[dz(t)dz*(t)]e- S , S = J dt(z(t)z*(t) - z(t)z*(t)), (A.15) 

where 

[dz{t)dz*(t)\ = [2dRez{t)dlmz{t)}. (A. 16) 



Appendix B: Another derivation of (14.81) 



In this appendix, we give another derivation of (|4.8|) . We first consider a Fourier transform 

of U(z, z*,t) 



U(A,A*,t) = J dzdz*e A * z - Az *U{z,z*,t), (B.l) 
which is an analogue of u. By substituting ()3.18|h we obtain 



Om+n ffm+n 

f/(A,A*,t)= Y ciC n e i{m - n)t \^(-l) m — — — e"2 AA *. (B.2) 



. m\n\ dA m dA* n 

m,n=0 



Some algebra gives 



[U(A,A*,t),U(A',A'*,t)) 



(&i - e* AA )U(A + A', A* + A ,t) 
e i ( A+A0(A* + A'*)-iAA*-iA'A'* 2 . s i n i (A* A 7 - AA'*)U{A + A', A* + A'*,t) (B.3) 
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This result urges us to define a new object 

U new (A,A*,t) = e 1 + AA *U(A 1 A*,t). (B.4) 

In fact, U new (A, A* ,t) satisfies the Woo algebra ()2.27j) with the identification A = (3 — ia. 
It is also easy to see that under this identification U new (A, A* ,t) satisfies the same equation 
of motion (|2.25|) and constraint (j2.26|) as u(a,/3, t) satisfies. Therefore U new would coincide 
with u: 

u(a, (3, t) = e^ AA *U{A, A*, t). (B.5) 



We also checked this by explicit calculation. The algebra (|B.3|) and the relation (jB.5|) were 
already derived in in the context of the integer quantum Hall effect. Thus we obtain a 
direct relation between U and u: 

dAdA* _ A . 2+A2 »_i AA » f dpaq a ( q+ j p ) + *£. ( q -i P ). 



U(z,z*,t)= / rrr^ e -A-*+A*-- a AA- / ^L e ~-(i+v)+-^p)u^ q ^). (b.6) 

J ^ J 27T 

Using ()B.6|) . we calculate §j(t) as 



d J f _i AA » f dpdq _a 



dA* J V J 2tt 

The last expression actually gives (|4.8jl . 



(B.7) 

A=A*=0 
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